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ABSTRACT 

We accurately determine a new Cepheid distance to MlOl (NGC 5457) using archival 
HST/ACS V and / time series photometry of two fields within the galaxy. We make a 
slight modification to the ISIS image subtraction package to obtain optimal differential 
light curves from HST data. We discovered 827 Cepheids with periods between 3 and 
80 days, the largest extragalactic sample of Cepheids observed with HST by a factor 
of 2. With this large Cepheid sample we find that the relative distance of MlOl from 
the LMC is A/i^j^Q = 10.63 =b 0.04 (random) it 0.07 (systematic) mag. If we use 
the geometrically determined maser distance to NGC 4258 as our distance anchor, the 
distance modulus of MlOl is /io = 29.04 it 0.05 (random) it 0.18 (systematic) mag or 
D = 6.4 lb 0.2 (random) it 0.5 (systematic) Mpc. The uncertainty is dominated by the 
maser distance estimate (±0.15 mag), which should improve over the next few years. 
We determine a steep metallicity dependence, 7, for our Cepheid sample through two 
methods, yielding 7 = —0.84 it 0.22 (random) it 0.07 (systematic) mag dex~^ and 
7 = — 0.721q22 (random) it 0.06 (systematic) mag dex~^. We see marginal evidence for 
variations in the Wesenheit P-L relation slope as a function of deprojected galactocentric 
radius. We also use the TRGB method to independently determine the distance modulus 
to MlOl of /io = 29.05 ± 0.06 (random) ± 0.12 (systematic) mag. 

Subject headings: Cepheids — distance scale — galaxies: individual (MlOl) 

1. Introduction 

Hubble Space Telescope (HST) observations of Cepheid variables have proven to be invaluable 
for e xtragalactic distance estimates outside th e local group, determining distances out to ~ 35 Mpc 
(e.g. iFreedman etaPboOll : IrIcss et al.ll2009al lbh. Type la supernovae and other secondary distance 



indicators used at greater distances rely on calibrations using extragalactic Cepheid distances. 
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Through these cahbrations, measurements of the Hubble constant (Hq) depend directly on the 
uncertainties in the Cepheid distance anchor. Utilizing the Wide Field an d Planetary Camera 2 
(WFPC2) aboard HST, th e HST Key Project IPreedman et alJll994l . boOlh and the Supernova la 
HST Calibration Program (jSaha et al.ll200ll : ISandage et al.ll2006l ) were both able to determine Hq 
with < 10% uncertainties. The two projects found estimates of Hq = 72 ± 8 km s~^ Mpc~^ and 
Hq = 62.3 lb 1.3 (random) ±5.0 (systematic) km s^^ Mpc^"*^, respectively, and the measurements 
were only consistent at the 2-sigma level. 

A dominant source of systematic error for both these measurements is their dependence on the 
distance to the Large Magellanic Cloud (LMC) for the absolute calibration of the Cepheid distances. 
The distance to the LMC is only kno wn to a precision of ~ 10%, with some estimates differing by 
±0.25 mag (e.g. iBenedict et al.ll2002l ). With such a large uncertainty in the distance, it is important 
to replace the LMC as our extragalactic distance anchor. There have been so me attempts to provid e 
a replacement anchor, for example the eclipsing binary dista nce to M33 of iBonanqs et al.l (j2006l ). 
NGC 4258, with a geometrically determi ned maser distance (iHerrnstein et al.lll999l ). is one of the 
most promising alternatives to the LMC. Macri et al. J 200^) us ed HST and the Advanced Camera 



Ford et al 



20031 ) observations of Cepheids in NGC 



for Surveys/Wide Field Camera (ACS/WFC; 
4258 to make the maser host galaxy a viable distance anchor for extragalactic Cepheid distance 
measurements. 

Extragalactic Cepheid distance measurements assume the Cepheids are at a common distance 
and apply the period-luminosity (P-L) relation of Cepheids variables. These measurements are 
complicated by the compositional dependence of the Cepheid P-L relation. The magnitude of this 



ory ( 


Stothers 1988: 


Stiftlll990: Fiorentino et al.ll2002: \ 


'^alle et al.ll2009 ) as well as Galactic and extra- 


galactic observations 


(Caldwell & Coulson 


1986. 


1987: 


Freedman & Madore 199C 


Gould 


1994: 


Stift 


1995 




Sasselov et al. 


1997: Kochanek 1997 


: Kennicuti 


19981: Groeneweecen et al. 


2004: ISakai et al. 


2004 


Macri et al. 


2006). Furthermore. Tammann et al. (2003) and Sandaee & Tammann ( 


2008) 



(and references therein) argued that th e slope and zero-point of the P-L relation must be dependent 
on metallicity, although iMadore et al.l (|2009l ) demonstrated that the slope of the "reddening free" 
Wese nheit P-L re l ation is insensitive to changes in the slope of the underlying V and / P-L rela- 



tions. 



Riess et al.l (l2009bl ) do not see any dependence in the slope of fit Wesenheit P-L relations in 
a sample of 7 galaxies observed with WFPC2 and ACS/WFC with average metallicities spanning 
0.2 dex. 

The purpose of this paper is to measure the distance to MlOl and test the universality of the 
Wesenheit P-L relation. We employ ISIS image subtraction (lAlard &: LuptonI 119981 : lAlardI 12000) 
on archival F555W and F814W time series ACS/WFC data of MlOl to identify a large sample of 
Cepheids. Using NGC 4258 as our distance anchor we are then able to measure the distance to 
MlOl. With our large sample of Cepheid variables we explore possible variations in the slope and 
zero-point of the Wesenheit P-L relation as a function of deprojected galactocentric radius. We 
investigate the possibility of a constant slope and zero-point, a linearly varying zero-point with a 
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constant slope, a linearly varying slope with a constant zero-point, and a linearly varying slope and 
zero-point. 

The paper is organized as follows: Section [2] describes the observations; Section [3] details the 
data reduction and image subtraction; Section H] explains our Cepheid search algorithms; Section [5] 
presents our Cepheid distance to MlOl; Section [6] investigates the universality of the Wesenheit P-L 
relation with deprojected galactocentric radius; Section [7] presents our Tip of the Red Giant Branch 
(TRGB) distance to MlOl; Section [8] shows interesting non-Cepheid variable objects; Section [9] 
discusses our results; and Section [10] presents our six main conclusions. 



Observations 



For this work we used archival ACS/WFC data (GO program 10918; PI Wendy Freedman) of 
MlOl (NGC 5457). MlOl has beer i classified as a SAB(rs)cd II I-IV spiral galaxy (jde Vaucouleurs et al 



1991 



) and a S(s)c I spiral galaxy (ISandaee fc Tammann 



1981 



The images were first processed by 
the STScI ACS calibration pipeline (see iPavlovsky et al.ll2005l ). We obtained the calibrated and 
flat-fielded images (_flt.fits) and the drizzled images (_drz.fits) from the Space Telescope Science 
Institute (STScI) HST Archive. The program observed two fields in MlOl with Field 1 centered 
at {a, 6) = (14'^03°'28!8, +54°21'35'.'0) and Field 2 centered at (14^02'^51!0, -F54°20'22'.'0). Fig. [U 
shows the position of these two fields as well as the WFPC2 fields used in the Hubble Key Project 
superimposed on the Sloan Digital Sky Survey (SDSS) image of MlOl obtained using the Hubble 
Legacy Archive^ 

The fields were observed in both F555W( F-band) and F814W(/-ba nd) for 12 epochs foL 



1994 



lowing a power law sampl ing in time to minimize possible period aliasing (jFreedman et al 
Madore fc FreedmanI |2005| ) . The observations have a baseline of ~ 30 days and were completed 



between 2006 Dec 23 and 2007 Jan 21 (Field 1) and between 2006 Dec 25 and 2007 Jan 23 (Field 
2). Table [T] provides a log of the observations. The observations were taken in a two point dither 
pattern with total exposure times at each epoch of 1330s for F555W and 724s for F814W. 



3. Data Reduction 

To make the most of these excellent data, we reduced them using a procedure that took 
advantage of both the increased sensitivity to variability that image subtraction provides and the 
high photometric precision inherent in instrument-dedicated photometry packages, such as the ACS 
module of the DOLPHOT package. Surprisingly, image subtraction has been underutilized on HST 
data even though it is well established in variable object searches at a variety of wavelengths (e.g. 



littp://hla.stsci.edu/ 
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Henderson et alj l20ld : iKhan et al 
2OO7I : 



Wozniak 



2010 



Koztowski et alJ I2OI0I : iNataf et alj I2OI0I : IPeeples et al 



2OOOI ) . Additionally, iBonanos Sz Stanekl (j2003l ) demonstrated the advantage of using 



image subtraction over more conventional photometric procedures in Cepheid variable observations, 
finding a nine-fold increase in the number of Cepheids found in Very Large Telescope observations 
of M83. To our knowledge, t his work is th e first paper applying image subtraction to HST data 
although the idea is not novel ( BersieilboOll ). The method we describe below uses image subtraction 
to produce light curves for all objects in our field in the instrumental magnitude system. We then 
use the ACS module o f the D OLPHOT package to measure and correct the magnitudes of these 
sources to the iLandoltl (| 19921 ) V and / magnitude system. 



3.1. Image Subtraction 



We used the image-subtraction package ISI^ ( Alard &: Lupt' 



ton 



1998 



Alard I2OOOI ). ISIS per- 



forms image subtraction utilizing a spatially variable kernel, allowing for accurate diff'erential pho- 
tometry even in crowded fields. 

Before we ran ISIS, we used standard IRAfU routines to multiply the drizzled images down- 
loaded from the HST archive by their exposure time to change their units from flux to counts 
so that we can accurately determine the uncertainties in the photometry. We then interpolated 
individual frames to a single drizzled image in each field, so all exposures had the same astromet- 
ric positions. This is usually done with ISIS, but ACS data is especially difficult as it has large 
numbers of resolved sources. ISIS was originally developed for microlensing and assumes that all 
sources are point sources. Instead we used Sexterp (Siverd 2011, in preparation)0, a hybrid code 
which utiliz es the interpolation met hod of ISIS with source identification and centroiding done by 
Sextractor (jBertin &: Arnoutslll996l ). Using Sexterp we interpolated all the images to match the 
first F814W image in each field (j9o401010_drz.fits for Field 1 and j9o413010_drz.fits for Field 2). 
As a result, all of our quoted astrometric coordinates will be referenced to the header astrometry 
of these two images. 

We used ISIS to create a reference image for each band in each field. First, ISIS transforms 
all the images to the same point-spread function (PSF) and background level by convolving the 
images with a space-varying convolution kernel. ISIS then stacks the resulting images using a 3- 
sigma rejection limit from the median. Fig.[2]and Fig. [3]show color composite images of the F555W 
and F814W reference images. 



We then used DAOPHOT and ALLSTAR (IStetsonlll987l . 



1992 



19941 ) for photometry on the 



•^The ISIS package and a tutorial are available at littp://www2.iap.fr/users/alard/package.litml 

■^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 

^http://www.astronomy.ohio-state.edu/~siverd/soft/is3/index.html 
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four ISIS reference images, to find more then 250, 000 objects on each image. We then matched 
objects in each field between the two filters and found that 142, 018 and 160, 551 matched in Field 
1 and Field 2 respectively. 

We next used ISIS to subtract each individual epoch from the corresponding reference image 
by first convolving the reference image with a space- varying kernel to match it to the image before 
subtracting it. In these subtracted images, all constant sources will disappear, leaving only the 
variable objects. ISIS makes a variability image by taking the mean absolute deviation of the 
subtracted images. A small region of the reference image and the variability image are shown in 
Fig.H 

We then used ISIS to extract the light-curves of all the sources identified in the reference images 
by DAOPHOT. However, the resolved sources inherent in HST observations cause ISIS to make a 
poor PSF itself. Instead we fed ISIS the PSF created in DAOPHOT using carefully selected isolated 
stars. This procedure requires a slightly altered version of ISIS that we describe in Appendix [Al 
Using this model PSF, we performed photometry on all the objects matched between the filters on 
all the subtracted images. ISIS outputs the light curves i n differential fiux unit s which we convert 



to an instrumental magnitude system using VARTOOLS ([Hartman et al.ll2008l ). 



3.2. DOLPHOT magnitude calibration 



The ACS module for DOLPHOT 1.10 was used on the calibrated and flat-fielded images of 
MlOl following to the recommended settings in the DOLPHOT User's GuidcH and DOLPHOT 
ACS User's Guid^. Photometry for each epoch and each field was performed separately with the 
same astrometric reference images as were used for ISIS. The recommended cuts were applied to 
remove non-stellar detections and to eliminate objects with uncertain photometry. For each epoch, 
we only accepted objects for which DOLPHOT reported a magnitude in both filters. After these 
cuts, DOLPHOT reported ^ 320, 000 objects in Field 1 and ^ 380, 000 objects in Field 2 for each 
epoch. 

The output of DOLPHOT contains an instrumental and transformed magnitude for each fil- 
ter as well as measurements for each of the individual, dithered frames. We used the combined 
transformed magnitudes for each filter. The final magnitudes ar e corrected for charg e transfer effi- 
ciency (CTE) losses. DOLPHOT uses transform ations given in Sirianni et al.l (j2005l ) to transform 



from the HST filters F555W and F814W to the iLandoltl (jl992l ) V and / magnitudes. However, 



Sirianni et al.l (|2005l ) warn that these transformations have a large uncertainty and that theoretical 



'''http://purcell.as.arizona.edu/dolphot/ 

®http: / /purcell. as. arizona.edu/dolphot / dolphot.ps. gz 

^http: / /purcell. as. arizona.edu/dolphot / dolphot ACS.ps.gz 
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and empirical transformations differ. Unfortunately, these transformations are required to make 
use of any published Cepheid P-L relations. Many systematic biases can be avoide d when using the 



maser-host galaxy NGC 4258 as a distance anchor, since the iMacri et al.l (|2006l ) Cepheid sample 
was also observed with HST ACS/WFC and applied the same transformations. In such comparisons 
these photometric issues cancel. 

To calibrate the light curves from image subtraction, we matched the DAOPHOT objects to 
DOLPHOT objects in each epoch, making sure to note the offset in positions due to the differing 
coordinate conventions between the two packages. We accepted DAOPHOT objects that matched 
an object in at least 1 DOLPHOT epoch, leaving 130,522 objects in Field 1 and 141,817 objects 
in Field 2. Differences between the DOLPHOT transformed magnitudes and ISIS instrumental 
magnitudes were then calculated individually for each object at each epoch where DOLPHOT 
had a matching object. We then used a weighted mean with iterative sigma clipping to find the 
magnitude offset for each object. Finally, each light curve was calibrated with the addition of 
this offset magnitude. Uncertainties in this offset were then added in quadrature to the Cepheid 
template fit mean magnitudes presented in §4.21 



Table [2] contains the positions and the calibrated V and / magnitudes of bright non- variable 
secondary standards. The magnitudes were measured on the first epoch for each field. 



4. Cepheid Variable Search 



In the following subse ctions we de s cribe our variability sea.rch an d Cepheid selection criteria, 
which is similar to that of Macri et al. ( 2006 ) and Riess et al. ( 2009bl ). A variability index is first 
used to drastically cut the number of sources. Light curves a re then fit wit h a Cep heid template 
and poor fits are determined by an F-test and rejected. The lUdalski et al.l (jl999al ) P-L relations 
were adopted and color cuts applied to refine the Cepheid sample. Lastly, a minimum period cut 
and an iterative distance modulus fit removed Population II Cepheids and other outliers. 



4.1. Variability Search 



The IWelch &: StetsonI (jl993l ) variability index /y was calculated on each of the light curves. 



The variability index is given by, 



\/ n[n — 1) 



where 5Vi and 6Ii are the residuals from the mean V and / magnitudes normalized by measurement 
uncertainties. Correlated deviations from the mean magnitude will contribute positively to the 
variability index whereas anti-correlated deviations contribute negatively. Thus, a constant star 
with random noise will have a variability index of ~ and a truly variable object will have a 
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positive variability index. Fig. [5] shows the variabihty index as a function of V magnitude for Field 
1. We then selected objects with ly > 1.5 as variable objects. Lastly, we removed objects with 
mean (V — I) < mag and mean {V — I) > 2.5 mag. These color cuts are well outside the expected 
range of the instability strip and lessen the number of variable objects that had to be fit with a 
Cepheid template. 



4.2. Cepheid Template Fit and F-Test 



curves 



We then fit the C epheid templates of lYoachim et al.l (|2009l ) to the variability-selected light 
from The lYoachim et"al] (I2OO9I1 templates are a principle coinponent analysis of a 



Four i er decomposit i on of Galatic Cepheids ( Ochsenbein et al. 



2001 



Gieren 



Henden 



1981 



1996 



Sebo et al 



Moffett fc Barnes 



Barnes et al 



2002 



Moffett et al 



1984 



2000; 



Coulson &: Caldwell 



Berdnikovlll997 : Berdnikoy Sz Turner 



1985 



Berdnikov Sz Turner 1995; 



199711. LMC Cepheids a n d SM C Cepheids (jUdalski et al.l Il999ai ib: 



1998). Yoachim et al. ( 2009 ) provide^ an Interactive Data Lan- 



guage (IDL) procedure to fit these Cepheid templ ates to light cu rves. As the amplitudes of different 
Cepheids with the same period are known to vary (IStetsonlll996l ). we modified the IDL Cepheid code 
to allow the amplitude of the Cepheid templates to vary by a multiplicative factor. Naively, one 
would allow the first principle component of the Cepheid template to be left as a free p arameter in 
the fit, as the majority of this component is amplitude variations ( Yoachim et al. 2OO9I1. However 



we fou nd that even when bounded to values within the scatter seen in Figure 8 of 
(|2009l ). this resulted in unrealistic Cepheid light curve fits. 



Yoachim et al 



The Cepheid fitting c ode implemented t he IDL Levenberg-Marquardt least-squares minimiza- 
tion procedure mpfit.prj^ ( Markwardt 2009 ) to fit the light curves giving an initial period and an 
initial phase. We started the fitting procedure at the two periods identified by the Minimum String 



Length Method (jBurke et al.lll97d : 



DworetskyI 1 19831 ) for the V and / light curves. However, no 



periodicity searching method will work for time baselines longer then the observational baseline. 
Fortunately, the accuracy of HST photometry and the use of Cepheid templates allows us to deter- 
mine the properties of Cepheid variables with periods significantly longer than this baseline. This 
requires a brute force template fitting search for periods longer than the baseline. In addition to the 
Minimum String Length periods, we also tried 250 initial periods from 3-100 days logarithmically 
spaced for each Cepheid. The lowest initial period bound of 3 days we set because Cepheids with 
shorter periods would be too close to our detection limit to be useful. We set an upper b ound of 100 
days i n the initial periods to avoid the possible change in the P-L relation as shown in bird et al] 
(|2009l ) for ultra-long period Cepheids. For each initial period we tried 10 different initial phases 
of the Cepheid template. We accepted the overall best-fit template for each variable based on the 
fit's x^value. Our method is computationally expensive but thorough. Representative light curves 



'http://www.astro.washington.edu/users/yoachim/Ceph_code.tar.gz 



'Available at |http : / / www. physics . wise . edu / ~craigm / idl/fitting . html | 
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for Cepheids in our final sample (as described in ^4.4p are shown in Fig. [6] and Fig. [7] for Field 1 
and Field 2 respectively. 

We then performed a version of the F-test to reject non-Cepheid variable objects. We compared 
the x^/dof of the Cepheid template fit (x^/dof template) with that of a constant magnitude least- 
squares fit (x^/dofeonstant) and a linear least-squares fit (x^/dofunear) for both V and / light curves. 
We then define, 

p _ min(x^/dofcontant,X^/dof linear) ^2) 
X^/dof template 

A stringent F > 3 cut was then applied. Fig. [S] shows the distribution of -F as a function of fit period 
for Field 1. The objects that survive this cut are the Cepheid candidates presented in Tables [3] and 
m along with their positions, variability index, fit period, fit mean magnitude and fit amplitudes. 



4.3. Fiducial Period-Luminosity Relations 



2001 



We adopted the updated OGLE-II P-L relations (jUdalski et al.l Il999al ) (as upd ated on the 
OGLE websit which have been the standard choice in most extragalactic studies (e.g. Freedman et al 



Macri et al.ll2006l ). Because we used the d i stance modulus to the LMC determined from the 



maser host NGC 4258 presented in lMacri et al.l (j2006l ). we adopted the same P-L relation for our 
final distance measurement. This choice preve nts possible systemat ic errors that could be caused 
by the adoption of a different P-L relation. The lUdalski et al.l (|l999al ) P-L relations were calculated 
on a sample of ~ 650 Cephieds in the LMC, and are given by: 



Vo = 14.287(0.021) - 2.779(0.031) [log P- 1] 
Iq = 13.615(0.014) - 2.979(0.021) [log P - 1] 



(3) 
(4) 



where P is the Cepheid period in days. We can then define the V and / LMC relative distance 
modulus: 

Anv = V-Vo, and Afij = I - Iq (5) 
where V and I are the measured Cepheid magnitudes. We then create the fiducial Cepheid color: 

{V - I)o = Vo-Io (6) 
using the adopted P-L relations. We can now define the color excess, 

E{V-I) = {V-I)-{V-I)o. (7) 



The LMC relative Wesenheit distance modulus (|Madordll982l ) is, 

A/XLMC=^W-(^-l)^(^-^)- 



(8) 



^ftp://sirius. astrouw.edu.pl/ogle/ogle2/var_stars/lmc/cep/catalog/README. PL. 
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where R = Ay /{Ay — Aj). Final ly, assuming a select ive to total extinction coefficient Ry = 3.1 



and using the Ax from Table 6 of lSchlegel et al.l (|l998l ) we calculate R = 2.45. 



4.4. Selection Criteria and Final Sample of Cepheid Variables 

We made the following cuts to our Cepheid candidate sample to identified true Cepheid vari- 
ables and to minimize the e ffects of blending and extinction. We chose these cuts following the 



criteria of lMacri et al.l (|2006l ) 



(1) Amplitude ratios - Objects with a ratio of / amplitude to V amplitude greater then 0.75 or 
less than 0.25 were removed from the sample. 

(2) Blue blends - Objects with E{V — I) < 2a below the Galactic foreground extinction value were 
removed from the sample. Cepheids that appear significantly bluer then the Galactic extinction are 
likely to have blue blends . We adopted E{B — V) = 0.009 mag or equivalently E{V — /) = 0.011 
mag ( Schlegel et al. 19981 ) as our foreground extinction. 



(3) Red blends and highly extinguished objects - Objects with E{B — V)> 0.5 mag or equivalently 
E{V — I) > 0.68 mag were removed from the sample. Cepheids that appear extremely red are either 
highly extincted or have red blends. Although the Wesenheit magnitude and Wesenheit distance 
modulus are resistant to reddening, highly extinguished sources may systematically differ if the 
true Rv value differs from our assumed Ry value. Additionally, if the color excess is due to blends 
these objects would clearly bias our distance. 

(4) Population II Cepheids - Objects with a distance modulus relative to the LMC of A/^lj^q > 
11.5 were removed from the sample to exclude W Virginis and RV Tauri variables. These Population 
II Cepheids have similar light curve shapes and colors as those of fundamental mode Cepheids, and 
some would have passed our initial sample selection. 

(5) Period Cut - Object with a period < 3 days were removed from the sample because their light 
curves are too noisy for accurate identification of Cepheid candidates. 

(6) Iterative Sigma Clipping - We then computed a weighted median LMC-relat ive distance 



modulus for each field using a least-absolute-deviation technique as also applied by iMacri et al 



(|2006l ). We perform iterative sigma clipping to remove outliers. This weighted median was used 
in place of the typical mean modulus because of its outlier resistance. Additionally, any Cepheid 
distance moduli distribution is skewed due to blending and the median will be a more robust 
statistic against these systematic problems. The difference between these two fitting methods is 
investigated in ^5.11 



The effects of the sample cuts are shown in Table [5l A period histogram for the final Cepheid 
sample is shown in Fig. [9l Cepheids and rejected Cepheids candidates are listed in Tables [3] and HI 
respectively. The quantities derived using our fiducial PL for the final Cepheid sample are listed 
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in Table [6j The locations of the Cepheids are shown in Fig. [10] and Fig. [TT] for Fields 1 and 
2, respectively. The color-magnitude diagrams (CMDs) of the fields with the Cepheids marked 
are presented in Fig. [121 The V, I and W luminosities for the final Cepheid sample are shown as 
functions of period with the predicted Wesenheit P-L relation overlain are shown in Fig. 1131 Finally, 
the V, I and W LMC-relative distance moduli are shown as a functions of period in Fig. [TH 



5. The Cepheid Distance to MlOl 



In the following subsections we describe how we determined the distance modulus of MlOl from 
our Cepheid sai nple. We first refined our Cepheid sample through cuts in period and color. We 
then assume the lBresolinI (j2007l ) oxygen abundance gradient for MlOl and compute the metallicities 
of our Cepheid sample. The distance modulus relative to the LMC was taken to be the distance 
modulus at [0/H] =8.50 dex of a linear least-squares fit to the ind i vidua l Cepheid distance moduli 
as a function of metallicity. Fi nally, we adopted the iMacri et al.l (j2006l ) LMC distance modulus. 



which was determined from the lHerrnstein et al.l (|1999 )~NGC 4258 maser distance, to compute our 
final distance modulus to MlOl. 



5.1. Refined Cepheid Sample for Distance Determination 

Because our goal is to determine the most accurate distance to MlOl, we require a refined 
sample of fundamental mode Cepheids. First, we exclude Cepheids with periods below 6 days 
because we are unable to distinguish between overtone and fundamental mode pulsators with such 
sparsely sampled light curves. Overtone pulsators with periods < 6 days have simi lar light curve 



shapes but are significantly brighter than their fundamental mode counterparts (e.g. lUdalski et al 



1999al ). and thus would bias our distance measurement. 



We computed the median distance modulus relative to the LMC under minimum period cuts 
from 6 to 20 days followin g the procedure des cribed in ^4.4[ This analysis is similar to the period 



cuts performed in §4.2.1 of lMacri et al.l (|2006l ) and we share the motivation therein. Fig. [T5] shows 
the median distance moduli relative to the LMC as a function of the minimum period for Fields 1 
and 2. We find that lower period cutoffs lead to smaller distance moduli. To avoid this bias, we 
adopt a minimum period of 14 days for both fields. Fig. [16] and Fig. [TTlshow the distance modulus 
relative to the LMC as a function of period for Fields 1 and 2, respectively. Fig. [18] and Fig. [19] 
show the E(V-I) color excess as a function of period for Field 1 and Field 2 respectively. These 
four plots also show the rejected Cepheid candidates from Table [4] for reference. These period cuts 
lead to a median distance modulus relative to the LMC of A^-^y^(j = 10.58ib0.01 mag using 155 
Cepheids in Field 1 and A^^j^^^ = 10.57±0.02 mag based on 135 Cepheids in Field 2. For each of 
the minimum period cuts, we also show the difference between the mean and median fit distance 
moduli relative to the LMC in Fig. [20] At a period cut of 14 days, there is a small difference in 
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the distance moduli O.OSitO.Ol mag and 0.03ib0.02 mag for Fields 1 and 2, respectively. 



For our final distance measurement C epheid samp l e we make more stringent cuts than were 
performed in ^4.41 as were also applied by iMacri et alj (j2006l ) to remove blue blends, remove red 
blends and ensure that our sample cover the same range in extinction. We discard Cepheids with 
color excesses below the Galactic foreground reddening E{B — V) < 0.011 mag, or, equivalently, 
E{V — I) < 0.015 mag. We also remove Cepheids that show large reddening, E{B — V) > 0.28 
mag, or, equivalently, E(y — I) > 0.39 mag, leaving 235 Cepheids in our refined Cepheid sample. 



5.2. Deprojection and Metallicity Dependence 

Extragalactic Cepheid studies make the implicit assumption that the galactic oxygen abun- 
dance gradient tracks the metallicity of Cepheids at the same galactocentric radius. With this 
assumption we obtain a differential measurement of the metallicity dependence of the Wesenheit 
P-L re lation. MlOl displays one o f the strongest radial abundance gradients found among nearby 
spiral ( Kennicutt fc G arnett^^ l996 ). We adopted the following oxygen abundance gradient for MlOl 
from Eq. 5 in iBresolin (i2007i ) : 



[0/H] = (8.75(0.05) - 0.90(0.07)p) dex 



(9) 



where p is th e fractional isophotal. iBresolinI ()2007l ) extended the 20 H ii oxygen abundance mea- 
surements of iKennicutt et al.l (120031 ) to within l'.5 o f the gala c tic ce nter with two additional mea- 
surements in the inner regions of MlOl. Figure 2 of Bresolin (j2007l ) shows the oxygen abundance 
gradient fit for MlOl where the H n regions span more then a dex in metallicity. We then com- 
puted the deprojected galactocentr ic distances f or the Cepheids were calculated using the following 
correcteJ^ equations presented in Macri et al. ( 2006 ): 



y 

P 



(a — ao) cos 5o cos {5 — Sq) sin 
{6 — 5o) cos (p — {a — ao) cos 5o sin ( 



cos I 



(x^ + 

Rq 



(10) 

(11) 

(12) 



where Rq is the isophotal radius. We adopted the nuclear position of MlOl from ICorwin et al 
(Il994l 'l of (J2000) {a, 6) = (14^03™ 12?53, +54° 20'55'.'6). We also adopted the maj or axis position 



m 



iior 

3). 



angle = 35?0 ± 1?0 and inclination of i = 17 ?0 ± 1?0 from IZaritskv et al.l (|1990l ). Additionally, 
w e assumed isoph otal radius Rq = 14'.4 from Ide Vaucouleurs et al.l (|l99ll ) as was also assumed 



BresolinI (|2007l ). Finally, we computed the metallicities of the Cepheids in our final sample as 



reported in Tabled 



Eq. 11 and Eq. 12 in Macri et al. 1 20061 ) are missing the geometric factor of cosiJo- 
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Fig.l21lshows the LMC relative distance modulus as a function of the deprojected galactocentric 
distance and metallicity for the refined Cepheid sample. A linear least-squares fit and a constant 
least-squares fit to the data are shown in Fig. [2TJ The standard deviation about the linear fit and 
the constant fit are 0.17 mag and 0.18 mag, respectively. The values of the linear fit and constant 
fit are x^=1077.1 and x^=1142.8, respectively. We see that the linear least-squares fit is preferred 
by the data. Fig. [22] displays the residuals from the linear least-squares fit against E{V — I), where 
no dramatic trends are seen. 

A linear least-squares fit to the data yields a slope of 7 = —0.84 it 0.22,. it 0.07^ mag 
dex~^ and a zero-point determination at [0/H] = 8.50 dex of A/x^j^q = 10.63 it 0.02^ ± 0.04^ 
maj^. We determined the random uncertainties of the linear best-fit parameters by bootstrap 
re-sampling the Cepheids with 2.5 x 10^ simulations. In each simulation we randomly re-sampled 
the refined Cepheid sample with replacement, to create a new sample with the same size. For 
each re-sample we calculated the linear least-squares best-fit slope and zero-point. We then took 
the standard deviation of the distribution of 2.5 x 10^ best-fit slopes/zero-points as the random 
uncertainty in the original determined slope/zero-point, respectively. As both fields in the data are 
at approximately the same deprojected galactocentric distance, we have a relatively small range in 
metallicity with which to determine the correlation between metallicity and luminosity, leading to 
the lar ge uncertainty w e see in 7. The systematic uncertainties are calculated from the uncertainties 



in the iBresolinI (j2007l ) oxygen abundance gradient. The random and systematic uncertainties in 
the zero-point are included in the distance error budget shown in Table [TJ represented by R5 and 
S4, respectively. 

To test if the Cepheids at the smallest galactocentric radii dominated the linear least-squares 
fit, we removed Cepheids with metallicities less than 8.65 dex from our sample and refit. The new 
best-fit slope 7 = —0.66 =b 0.25^ it 0.05^ mag dex~^ is shallower but consistent with the slope 
determined before the metallicity cut. Even with a large uncertainty in 7 there is tension (2.2 a) 
between our measured value and the Hubble Key Pro ject determination of 7 = —0.24 ±0.16 mag 



dex"^ for MlOl using WFPC2 data (lKennicuttill998l ). The implications of our measured Cepheid 
metallicity dependence are discussed in ^ 



5.3. Final Cepheid Distance to MlOl 

We then adopted the iMacri et al.1 kood) LMC di stance modulus ^LMC = l^-^l ± 0.10^ ± 
0.13s mag determined from the iHerrnstein et al.l (|l999l ) NGC 4258 maser distance to compute our 
final distance modulus to MlOl. Our error budget is shown in Table [71 We report a distance 
modulus relative to the LMC of A/^lj^jq = 10.63 ± 0.04^ =b 0.07s mag and a MlOl distance 



^ In this paper we denote random and systematic uncertainties with subscript r and s respectively, i.e., 3.00 ± 
O.lOr ±0.10. mag. 
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modulus of ^0 = 29.04 it 0.05^ ± 0.18s mag. We compare these distances to other determinations 
of distance modulus relative to the LMC and estimates of the distance to MlOl in Tabled 

6. Slope and Zero-Point of the Wesenheit P-L Relation 

We can also use the large number of extragalactic Cepheids in our sample to test the univer- 
sality of the Wesenheit P-L relation. Let us define the classic "Wesenheit reddening-free" mean 
magnitude by: 

W = V - R{V - I), (13) 

where R is defined in ^4.31 Based on the mean V and I magnitudes of our final Cepheid sample 
given in Table [6] we can define a expected mean Wesenheit magnitude (Wq) for each Cepheid. We 
then want to fit for a Wesenheit P-L relation: 

Wo = a^logP + b^ (14) 

where the slope and zero-point may depend on deprojected radius as a proxy for metallicity. 
Based on Eq. [HI the expected mean Wesenheit magnitude as a function of period and p for each 
these models is given by: 

WoiP) = [al + aiip - p')]ilogP- logP') + [hi + b^ip - p')] (15) 

where and 6^ are the values of the slope and zero-point, respectively, of the Wesenheit P-L 
relation at some characteristic fraction of the isophotal radius (p'). and b^ are the rates at 
which the slope and zero-point of the Wesenheit P-L relation change with p. The zero-point of the 
Wesenheit P-L relation is defined at a characteristic pulsation period (P')- order to minimize 
the covariance between a^, a^, 6^ and b^ we use the mean log period (logP' = 1.3 log(days) 
~ (logP)) and mean factional isophotal radius {p' = 0.2 ~ (p)) for the characteristic period and 
radius, respectively. 

We use the Cepheid sample in Table [3] and imposed the same minimum period cut of P > 14 
days used in ^5.11 and either impose no maximum period cut or impose a maximum period cut of 
P < 40.0 days to define a full and restricted sample, respectively. There are two reasons a maximum 
period may avoid systematic deviations. First, the ~ 30 day baseline of our data makes the period 
and magnitude determination suspect for very long period Cepheids. Many of them are missed 
in any case because they fail our F-test in §4.21 Eliminating the long period Cepheids is of little 
concern when determining the distance modulus to MlOl because the mean magnitude and period 
errors will also be larger on these objects. However, our ability to estimate the radius-dependent 
slopes of the Wesenheit P-L relation could be strongly affected by the stochastic nature that the 
inclusion of a very few long period Cepheids would induce. However, adopting this maximum cut 
limits our sensitivity to changes in the P-L slope. As there is no obvious solution to address both 
concerns, we performed and present both fits. After trimming, there are 297 Cepheids in our full 
sample and 278 Cepheids in our restricted sample. 
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We then fit Eq. [15] to both the fuh and restricted samples using the fohowing procedure. First, 
we fit a model with no radial dependence, and iteratively remove outliers. We took the 100 Cepheids 
closest in radius to each Cepheid, computed their scatter about the best-fit Wesenheit P-L, and 
rejected the Cepheid if it was more then 3 times the scatter from the best-fit P-L relation. This 
procedure takes into account the increased scatter of Cepheids at smaller p. We then renormalized 
the errors to have x'^/dof = 1 for the model with no radial dependence to account for the intrinsic 
scatter of the Cepheids. We then used the IDL Levenberg-Marquardt minimization procedure 



mpfit.pro (jMarkwardtl l2009l ) to fit four model Wesenheit P-L relations. Model 1 had a constant 
slope and zero-point. Model 2 had a linearly varying zero-point with a constant slope. Model 3 had 
a linearly varying slope with a constant zero-point and Model 4 had a linearly varying slope and 
zero-point. For each model initial and 6^ values were taken from the best-fit model with no 
radial dependence and the initial slopes and 6^ were taken to be 0. and 6^ were fixed to for 
the constant slope models and the constant zero-point models, respectively. To determine the 68% 
confidence intervals for each model we performed a bootstrap re-sampling with 2.5 x 10^ simulations. 
Fig. [23] and Fig. [26] present the Model 1, Model 2 and Model 4 fits on top of the Wesenheit P-L 
relations for six sub-samples with equal number of Cepheids stepping out in deprojected radius 
for the full and trimmed sample, respectively. The best-fit slope and best-fit zero-point, both as a 
function of p, for the full sample are shown in Figures [M] and [25] respectively. The best fit slope 
and intercepts of the six sub-samples in Fig. [23] are also shown for comparison in Figures [23] and 
respec tively. Best-fit parameters and values for the four model are presented in Table [9] 



Using the iBresolinI ([20071 ) abundance gradient we transform these fit parameter into functions of 
metallicity as shown in Table [10] The characteristic radius is then equivalent to a characteristic 
metallicity of 8.57 dex. 

We then compared the models and tested the validity of the additional terms between models 
fit to both samples. Based upon the x^values presented in Table [9] the data seems to prefer models 
2 and 4 over model 1 for both samples. We then computed the F-test confidence to quantify our 
justification for additional parameters between models 1 and 2 and between models 2 and 4. The 
calculated confidence levels are shown in Table [TT] Both Cepheid samples are better fit by Model 
2 than Model 1. This is expected because a changing Wesenheit P-L zero-point (6^) is equivalent 
to the change in the distance modulus as a function of metallicity (7) presented in Fig. [2T] When 
comparing Model 2 and Model 4 we note a marginal detection of a change in the Wesenheit P-L 
relationship slope with p. This detection is at a 93.6% confidence level for the full sample, but the 
significance decreases to 89.4% for the restricted sample. We postpone further discussion until §[9] 



TRGB Distance to MlOl 



We also determined the distance to MlOl using the TRGB meth od (iLee et al 



underlying physics of this method is well under stood and described i n 



effects of metallicity are also well understood ([Bellazzini et al. 



2001 



S alar is et al 



The 
The 

20041 ) ■ and many methods 



1993). 
(I2OO2I ). 
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exist to correct for these effects (jMager et al.ll2008l : iMadore et al.ll2009l ). The TRGB method must 
be apphed at the largest possible galactocentri c radii to liinit con tamination from the intermediate 
mass asymptotic giant branch (AGE) stars (jSakai et al.l 119961 ) . AGB stars have a continuous 
distribution at similar colors but extend to higher luminosities than the TRGB, diluting the abrupt 
edge of the old, low mass stars at the TRGB. 

The TRGB method works best given deep and accurate photometry. To accomplish this we 
reran DOLPHOT for each field, but include all epochs in a single run. DOLPHOT then combines 
each epoch's photometry and computes one transformed magnitude per object per band. This 
second run of DOLPHOT yielded 3.75 x 10^ and 3.90 x 10^ sources in Fields 1 and 2, respectively. 
We then trimmed the sample for optimally application of our TRGB detection method. We first 
eliminated sources at galactocentric radii smaller then 4.75' to limit the contamination from AGB 
stars, yielding 5.64 x 10^ and 3.69 x 10^ sources from Fields 1 and 2, respectively. We removed 
sources with y — / < 1.00 mag as they are too blue to be TRGB stars. Our final TRGB sample 
size contained 2.65 x 10^ and 2.22 x 10^ sources from Fields 1 and 2, respectively. 

We then appli ed reddening c o rrecti ons to account for the small amount of foreground extinc- 
tion. Based on the ISchlegel et al.l (| 19981 ) dust maps, the foreground absolute extinction correction 
of V and / are Ay = 0.028 mag and Aj = 0.017 mag. 



We then used these corrected magnitudes to calculate the lMadore et al.l (|2009l ) T magnitude 
for all our TRGB sources: 



r = /o-/3[(y-/)o-7TRGB], 



(16) 



where f3 is the / versus V — I slope of the TRGB, and 7trgb is the fiducial color where the T 
magnitudes are normalized to the / magnitude of a TRGB star of that color. The T magnitude 
is constructed to correct for the slope of the tip / magnitudes as a function of ^ — /, forcing 
the tip T magnitudes to be constant over a large range of colors or equivalently metallicitiea^^l. 
With the TRGB slant rem oved, edge de t ection methods more efficiently determine the TRGB 
position with fewer sources. Madore et al. (j2009l ) adopt (5 = 0.20 ±0.05 for the intrinsic color range 
1.5 < {V — I)o < 3.0, which is the color range o bserved for our T RGB. This slope is consistent 
with both observational valu es of /3 = 0.22 it 0.02 (jRizzi et al.l 120071 ) and theoretical predictions of 



= .15 (iMager et al.ll2008l ;i. With a fiducial color of 7trgb — ~ I)o — 1.5 mag, I Madore et al 



([20091) define the distance modulus based on the TRGB by: 



/J-TRGB = T — Mr 



T + 4.05(0.12) mag, 



(17) 



where Mtrgb is the TRGB absolute / magnitude at the fiduci al color 7. The uncer tainty in M^rgb 



was taken from Maeer et al. ( 20081) . a, com pani on paper to Madore et al. ( 2009 ). The absolute 



TRGB adopted in both lMadore et all (|2009l ) and iMager et al J 



same fiducial color, so we adopt the uncertainty presented in lMager et al.l (j2008l ). 



20081') are . 01 m ag different at the 



See Figure 1 in lMadore et al. for a schematic representation 
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The TRGB magnitude is determined using a edge detection procedure similar to that of 



Sakai et alj (j2004l ). but using the newly constructed T magnitude in place of the typical / mag- 
nitude. A c ontinuous luminos ity function (f)[m) was first com puted using equation (Al) in the 
appendix of ISakai et al.1 (|l996l ). We then used equation (4) of iMendez et alJ (|2002l ) to compute 
a logarithmic edge-detection response function E{m). This function is a continuous version of a 
Sobel filter. Similar to what is done in discrete edge detectors, w e weight the edge-det ector response 
by the Poisson noise of our luminosity function, E{m)[(j){m)]2 (IMendez et al.ll2002l ) and took the 
maximum as the TRGB. We then estimated the random uncertainty of tip determination by boot- 
strap re-sampling of the TRGB sample with 2.5 x 10^ simulations. In each simulation we randomly 
re-sampled the TRGB sample with replacement to create new sample of the same size. For each 
re-sample we determine the TRGB T magnitude with our edge detection method. We then took the 
standard deviation of the distribution of 2.5 x 10^ TRGB T magnitudes as the random uncertainty 
in the original tip detection. Fig. [27] presents the T magnitude CMD in the left panel with the 
tip detection response function in the right panel. The tip detection at T = 25.00ib 0.06 mag is 
marked on both panels for reference. With this detection, we use Eq. [T7]to determine the TRGB 
distance modulus of /xq =29.05±0.06r=t0.12s mag. The systematic uncertainty is dominated by the 
uncertainty of Mtrgb • 



8. Non-Cepheid Variable Objects 

There have been many studies searching for luminous variable stars, including a recent study 
of the LMC which f ound 117 variable objects with absolute magnitudes of My > —4.4 mag 



([Szczygiel et al.l 120101 ). We performed a non-exhaustive (by eye) search of the brightest, most 
variable objects in both fields which are not in our final Cepheid sample shown in Table [3j The 
majority of these objects appeared to be rejected Cepheid candidates, which usually displayed an 
obvious reason for their rejection, namely blends from colors or luminosity. Additionally, many 
possible long period Cepheids were rejected since they had poor light curve coverage due to the 
restricted time span of the data. 

Some of the most luminous and interesting variable objects are shown in Fig. [28] with their 
properties given in Table [T2l Possible periodic objects are shown in Fig. [29] with their properties 
listed in Table [T3l Fig. [30] presents the CMD of both fields with the luminous variables marked. 



9. Discussion 

In ^ and S}7] we determined the distance to MlOl using Cepheids and the TRGB method, 
respectively, and found that they are in good agreement with //coph = 29.04 it 0.05,. ± 0.18^, mag 
and /iTRGB =29.05±0.06r.±0.12s mag. We compared our MlOl distance moduli against other recent 
studies in Table [8] We see that our LMC relative distance modulus is in excellent agreement with 



-17- 



the recent measurements of I S aha et al.l (|2006l ) and iFreedman et alj (|200ll ) and our absolute MlOl 
distance moduli differ within the uncertainties. Thi s difference is due to the applied luminosity 
calibration for the adopted Cepheid P-L relations. ISaha et al.l (|2006l ) used a hybrid method to 



determine the ni e tallici ty dependence of the P-L r elation determi n ed us ing Galactic Cepheids from 
Tammann et al.l (|2003l ) and LMC Cepheids from ISandage et al.l (12004 ) . However, the Wesenheit 



P-L slopes measured for the sev en galaxy sample of 



than the Wesenheit P-L slope of lTammann et al 
difficulties with the method of the 



(20O 



Riess et al 



Tammann et al 



Riess et al. 



(l2009bl^ are significantlv shallower 



■_ .20091^ ) points out the many 

Galactic Cepheid P-L determination 



including the use of non-parallax based distances and high mean extinctions. 

Despite our large sample of Cepheids, our MlOl distance modulus measurement suffers from 
lar ge systematic unc e rtaint ies. As seen in Table [7l the majority of this uncertainty originates from 
the lHerrnstein et al.l (119991 ) maser distance to NGC 4258. Efforts to improve the uncertainty of the 



maser distance to NGC 4258 to 3% are u nderway using Very Long Bas eline Array obser v ation s 
(jArgon et al.l 120071 : [Humphreys et al.l 120081 ). However, to our knowledge, iHerrnstein et al.l (119991 ) 
is the latest actual measurement in the literature. Table [7] contains a goal error budget calculated 
under the assumption of a 3% maser distance to NGC 4258. Additionally, as this is a differential 
meas urement, systemati c errors in the photometric calibrations could be negated with a reanalysis 
of the lMacri et al.l (|2006l ) Cepheid observations of NGC 4258 with the same photometric procedure. 



We note that the Hq measurement with a total uncertai nty of 4.8% from lRiess et al.l (|2009al ) 
uses the published maser distance of /Xn4258 = 29.29 mag from IHerrnstein et al.l (jl999l ) but assumes 
the error bars (3%) of the yet unpublished result of Humphreys et al. (2009) rather than the 
published uncertainties (7%) of IHerrnstein et al.l (|l999l ). 



The universality of the Wesenheit P-L relation has been a point of concern and contention 
in the literature since it has important implications for deriving extragalactic distances. With 
our large sample of MlOl Cepheids we investigate the Wesenheit P-L as a function of deprojected 
galactocentric radius, and by assumption, as a function of metallicity. We determine the metallicity 
dependence of the Wesenheit P-L zero-point in §5.21 and Jj6] to be 7 = —0.84 it 0.22,. ± 0.07^ 
mag dex~^ and 7 = 6^ = — 0.72^0 22 (^^-iidom) ±0.06^ mag dex~^, respectively. These metallicity 
dependences are in contention with current values in the literature. Using the TRGB method as 
a metallicity independent distance indicator ISakai et al.l (j2004l ) report 7 = —0.25 ± 0.09 mag. 



However, the lSakai et al.l (|2004l ) results have recently be en challenged in iRizzi et al.l (j2007l ). The 
HST Key Project reported 7 = -0.24 ± 0.16 mag dex"! (|Kennicuttlll998l ^ using an WFPC2 data 
for an inner field and an outer field in MlOl. There is even contention with the more recent 

-0.29 ± 0.09,. ±0.05^ mag dex'^ using HST 



determination from lMacri et al.l (12(30611 who found 7 
ACS data of NGC 4258. However, iBono et al.l (j2008l ) shows that by removing a single outlier H 11 
region used in the determination of the metallicity gradients of NGC 4258 leads to a far shallower 
metallicity gradient: 



[0/H] = (8.89 - 0.16(p - 0.4) dex. 



(18) 
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then assumed bv lMacri et alj (|2006l ): 



[0/H] 



.97 ± 0.06) - (0.49 ± 0.08) - 0.4) dex. 



(19) 



With the shallower metallicity gradient of iBono et al.l (j2008l ). the 7 measurement of iMacri et al 
( 20od ) for NGC 4258 steepens to 7 = —0.89 mag dex~^, and is consistent with our estimate for 
MlOl. 

There have been arguments that the change in the zero-point of the Wesenheit P-L as a 
function of deprojected galac tocentric ra d ius ni ay be partially due to changes in blending as well 
as metallicity. For example, I Bono et al.l (|2008l ) argues that blending could explain the negative 
values of 7 that are seen because Cepheids in the inner region will be more blended on average, 
decreasing their apparent distance modulus. The idea that blending could undermine the commu- 
nity's efforts to use Cep h eids as a distance indica tor h as been disc u ssed i n many studies including 
Stanek fc Udalskil (|l999l ). iMocheiska et al.1 kood ) and lMacri etaP (l200lh . 



In addition to the varying zero-point of the We s enhei t P-L, the slope is also expected to 
show variations at some level. ISandage &: TammannI ([200a) argue that slope of the F-band P- 
L relation decreases monotonically with metallicity. However, iMadore et al.l (|2009l ) show that the 
Wesenheit P-L slope is largely insensitive to changes in the slopes of the underlying monochromatic 
P-L relations. Throug h an algebraic argum ent and reasonable assumptions about the observed 
properties of Cepheids, iMadore et al.l (|2009l ) show that the change in slope of the Wesenh e it P-L 
could be ~ — 0.3 times the change in the slope of the V P-L. However, IMadore et al.l (|2009l ) 
notes that the observational properties used are not currently known to a 10% precision, so their 
conclusions must be taken with caution. 

Our Cepheid sample is marginally better fit by a Wesenheit P-L relation model with both 
slope and zero-point linearly varying as a function of deprojected isophotal radius than by a model 
with a linearly varying zero-point alone. The evidence for a linearly varying slope as a function of 
p degrades when using a sample with an applied maximum period cut of 40.0 days. However, this 
degradation was expected as the shortest and longest period Cepheids are those that best constrain 
the slope of the Wesenheit P-L. We see three possible explanations for our observations: 

1. If Cepheid blending is a larger effect at smaller radii p, then one would naively expect a 
change in slope with our observed sign. To illustrate this point, imagine adding a constant flux to 
all Cepheids as a function of p; one would then expect more than just a change in the zero-point of 
the Wesenheit P-L relation. At smaller p the less luminous, shorter period Cepheids would increase 
in magnitude more then their more luminous, longer period counterparts at the same p. Thus 
the slope of the Wesenheit P-L relation would become shallower in the inner regions of a galaxy. 
This type of blending is not trivial or obvious, as the smooth increase in the stellar density would 
simply be subtracted in our photometric procedures. An increase in stellar density associated with 
Cepheids as a function of p would have to be invoked. A detailed and quantitative analysis of the 
magnitude of this effect is beyond the scope of this study. 
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2. From the recent work of ISandage &: Tammanij (|2008l ) and iMadore et al.l (|2009l ). it is ex- 
pected that the metalHcity dependence of the Wesenheit P-L slope should be small. First, assume 
that the change in the slope of the Wesenheit P-L relation is a negative fraction of the change in 
the slope of the V P-L relation (iMadore et al. 2009). Also a ssume that the V P-L slope decreases 
monotonically with metallicity (jSandage &: TammannI |2008| ) . From this we would expect the ob- 
served Wesenheit P-L relation to become shallower towards higher metallicities in the inner regions 
of a galaxy. 

3. The observed change in slope of the Wesenheit P-L relation is tenuous, and as such could 
simply be a statistical fluctuation in the data. 

MlOl is a large face-on spiral galaxy, has bright H ii regions at all radii and displays one o f 
the strongest radial abundance gradients found among nearby spiral (iKennicutt &: Garnettill996l ). 
which makes it an extremely useful object to study the systematic effects on the Wesenheit P-L 
relation for Cepheids. The two fields observed in MlOl used in this study are at approximately 
the same deprojected galactocentric radius, so the range in p is relatively small. This makes any 
study of the effects that p or metallicity have on the observed Wesenheit P-L relation difficult. It 
would be worth while to obtain another equivalent set of Cepheid observation of MlOl at a greater 
deprojected galactocentric radius to further test the effects of metallicity on the Wesenheit P-L 
relation. Additionally, it would be useful and of relatively little cost to reobserve one of the fields 
in MlOl used in this work for 2-3 additional epochs with HST/ACS. These observations would help 
to identify and constrain the periods of Cepheids with periods greater then 30 days. With a large 
sample of Cepheids covering a large range in period it would be possible t o test the chang e in the 
P-L relation slope for Ultra Long Period Cepheids at periods > 100 days (jBird et al.ll2009l ). 



10. Summary 

We accurately determine a Cepheid distance to MlOl (NGC 5457) using archival HST ACS 
V and / time series photometry of two fields within the galaxy. We used ISIS image subtraction 
to obtain differential light curves and developed a new method where ISIS uses an externally 
determined PSF model. We used DAOPHOT to create this PSF model and DOLPHOT to calibrate 
the photometry. Our six primary results are: 

1. We discovered 827 Cepheids using 12 epochs of HST ACS/WFC F555W and F814W data 
for two fields within MlOl. 

2. We determine a Cepheid based LMC-relative distance modulus to MlOl of A/iL]y[Q = 
10.63± 0.04,. ± 0.07s mag if we assume the OGLE II LMC F-band and /-band PL relations. 



3. Base d on the iHerrnstein et al.l (|l999l ) geometrically determined maser distance to NGC 
4258 and the iMacri et al.l (|20od ) Cepheid based LMC-relative distance modulus to NGC 4258, we 
determine an absolute MlOl distance modulus of po = 29.04 ± 0.05,.± 0.18^ mag. 
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4. We also determine a MlOl distance modulus based on the TRGB method of fiQ = 29.05 it 
0.06,. ± 0.12^ mag. 

5. We find a steep Cepheid metallicity dependence through 2 methods, yielding 7 = —0.84 it 
0.22^ ± 0.07^ mag dex'^ and 7 = 6^= -0.72^°;^° (random) ±0.06^ mag dex'^ 

6. Our Cepheid sample is marginally better fit by a P-L relation where both the slope and 
the zero-point linearly vary with galactocentric radius, over models with only a linearly varying 
zero-point. 

We would like to thank Christopher Kochanek, Lucas Macri, and Jennifer van Saders for 
their invaluable suggestions and discussions. We would also like to thank Joe Antognini, Thomas 
Beatty, Andrew Dolphin and Jill Gerke for other useful suggestions. This research gratefully made 
use of the NASA/IPAC Extragalactic Database (NED). Benjamin J. Shappee was supported by a 
Graduate Research Fellowship from the National Science Foundation. 



A. Creating a PSF for ISIS from Daophot 

Here we describe how to make ISIS use an externally determined PSF model. 

Before creating the ISIS PSF from a DAOPHOT PSF we must first use IRAF to change the 
units and zero-point of the ACS/WFC drizzled images. We must have run the modified ISIS 
interp.csh0, the stock script ref.csh and the stock script subtract. csh to create the ISIS reference 
images and subtracted images0 Additionally, DAOPHOT must be run on each reference image, 
carefully selecting a sample of PSF stars. Besides iteratively rejecting stars which differ from 
the DAOPHOT model PSF, we must also reject stars which have extended emission or objects 
with other stars nearby. Rejecting stars which have other nearby stars can be done using the 
WATCH PROGRESS option in DAOPHOT. We now have a PSF from DAOPHOT which is free of 
contamination from non-point sources which we save as PSF. fits. We also use standard procedures 
in DAOPHOT and ALLSTAR to create a catalog of source positions and photometry. We remove 
non-stellar sources from the catalog using cuts on the catalog. We then use the source list (.als) 
from ALLSTAR as an object list for ISIS by using the following awk command and removing the 
first 3 lines: 

awk '{printf "%9.3f 7.9. 3f 7.5. Of 7o5.0f 7ol3s 7.8. 3f 76 . 3f \n" , $2-1 , $3-1 , $2-1 , 
$3-1, "lc"$l". data", $4-0.1, $5}' ref.als > ! phot. data 



^*http: //www. astronomy.ohio-state.edu/~siverd/soft/is3/index. html 

^^The script detect. csh can also be run to create the stacked subtracted images seen in Fig. ID However, these 
images are not used to obtain the differential flux light curves in our analysis. 
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where the —1 pixel shift in the awk command accounts for the difference in the coordinate con- 
ventions between DAOPHOT and ISld^. This object hst should be the file associated with the 
VARIABLES parameter in the ISIS process_config file. 

We then need to manuahy make two files which ISIS usually creates and put them in the 
directory with the images to be processed by ISIS. The first is an nxn fits image named 'psfJileO. fits' 
filled with zeros, where n is specified by psf_width in the ISIS phot_config file. The second file must 
be named 'psf_table' which will contain one line of the form: 

psf_file0.fits 9 9 4211 4234 0, 

where the 1®^ element is the PSF file name, the 2^^^ and 3^"^ are the center of the 'psfJileO. fits' 
image in ISIS coordinates and the 4 and 5 are the x,y size of the reference images in pixels. 
We then use the addstar command in DAOPHOT to put a 'fake' star with extremely high signal 
to noise to the center of psf_file0.fits. It is important to account difference in the pixel coordinate 
system convention of DAOPHOT when creating the PSF. 

Lastly, the stock phot.csh script must be altered so that it does not call the program Bphot, 
which would usually make the ISIS PSF. The modified phot.csh script can then be run, creating 
the light curves for every object in the object list created in DAOPHOT. These light curves will be 
in differential flux units, but can be transformed to an instrumental magnitude using 'fluxtomag' 
in the VARTOOL^^ package. 
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Table 1. Log of Observations 
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Table 2. Secondary Standards 
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Note. — R.A. and Dec. are in units of degrees. Errors in mean magnitudes and somiamplitudes 
are in units of 0.001 mag. Table discussed in i|3.2l 



- 28 - 



Table 3. Cepheid Variables 
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22, 


.849(014) 


0.477 


0.278 


275.8 


Fl- 


■1767 


21.7926(0.3649) 


210.86424080 


54.34258320 


3383.2 


2107.0 


23, 


.669(021) 


22, 


.922(024) 


0.335 


0.202 


129.5 


Fl- 


■1890 


30.4055(0.7575) 


210.83646784 


54.36458112 


2380.2 


3798.2 


23, 


.699(015) 


22, 


.755(013) 


0.479 


0.269 


501.1 


Fl- 


■1975 


30.4193(0.5237) 


210.86646959 


54.34309782 


3312.5 


2035.4 


23, 


.728(016) 


22, 


.797(012) 


0.501 


0.266 


241.0 


Fl- 


■2045 


22.8373(0.1005) 


210.87188971 


54.34614160 


3021.9 


1911.6 


23, 


.776(013) 


23, 


.083(015) 


0.446 


0.302 


236.3 



Table to be published in its entirety in machine-readable form. 



Note. — R.A. and Dec. are in units of degrees. Astrometry presented is referenced to the HST headers of j9o401010_drz.fits for Field 
1 and j9o413010_drz.fits for Field 2. Errors in mean magnitudes and semiamplitudes are in units of 0.001 mag. Tabic discussed in i|4.2l 



Table 4. Rejected Cepheid Candidates 







Period 


R.A. 


Dec. 


X 


Y 




Mag 


nitudes 


Amp. 






ID 


(days) 


(J2000) 


(J2000) 


pixels 


pixels 




V 




/ 


V 


/ 


W 


Reason"' 


Fl- 


■381 


146.449 


210.86070136 


54.35415151 


2674.5 


2569.1 


22, 


,525(007) 


22, 


.265(011) 


0.136 


0.140 


98.5 


A 


Fl- 


■757 


166.814 


210.84322037 


54.38390199 


988.9 


4080.5 


22, 


.968(010) 


22, 


.648(014) 


0.281 


0.226 


62.6 


A 


Fl- 


■1005 


76.313 


210.87279394 


54.36153133 


1987.2 


2309.6 


22, 


,883(015) 


21, 


.845(010) 


0.336 


0.260 


422.4 


A 


Fl- 


■1357 


5.986 


210.86954527 


54.33204480 


3994.6 


1605.6 


23, 


.518(009) 


23, 


.128(014) 


0.120 


0.174 


40.6 


A 


Fl- 


■1378 


168.978 


210.88086502 


54.33618592 


3534.5 


1284.7 


22, 


.867(057) 


22, 


.828(070) 


1.279 


0.295 


19.2 


A 


Fl- 


■1576 


148.915 


210.88486096 


54.35099847 


2487.3 


1547.2 


23, 


.457(008) 


23, 


.252(016) 


0.139 


0.145 


21.0 


A 


Fl- 


■2754 


13.727 


210.88137081 


54.37066234 


1241.5 


2235.4 


24, 


.022(013) 


23, 


.661(019) 


0.249 


0.212 


62.8 


A 


Fl- 


-2812 


34.901 


210.86097991 


54.37769971 


1109.2 


3220.5 


23, 


.966(013) 


22, 


.849(009) 


0.314 


0.278 


135.8 


A 


Fl- 


■2883 


36.574 


210.87213007 


54.36013865 


2090.4 


2296.1 


23, 


,930(011) 


23, 


.003(012) 


0.374 


0.303 


236.1 


A 


Fl- 


■4020 


19.565 


210.85013743 


54.37919292 


1187.8 


3681.1 


24, 


.339(016) 


23, 


.475(017) 


0.190 


0.152 


25.2 


A 



Table to be published in its entirety in machine-readable form. 



Note. — R.A. and Dec. are in units of degrees. Astrometry presented is referenced to the HST headers of j9o401010_drz.fits for Field 1 
and j9o413010_drz.fits for Field 2. Errors in mean magnitudes are in units of 0.001 mag. Table discussed in ^4.21 

" Reasons for the Cepheid canidates rejection; A: Amplitude ratio ratio cut. B: Blue blends. C: Large Extinction. D: Population II 
Cepheids. E: Candidates removed by period cut. F: Candidates removed by iterative sigma clipping. 
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Table 5. Sample Selection 





Field 




Sample Cuts 


Fl 


F2 


Initial Cepheid Candidates .... 


632 


595 


1. Amplitude Ratio 


-139 - 


'121 


2. E{V ~ I) < (2o- foreground) 


-29 


-19 


3. eIv -I) > 0.5 mag 


-24 


-11 


4. Afj,o > 11.5 mag 


-6 


-7 


5. Period < 3 days 


-2 


-4 


6. cr clipping 


-21 


-17 


Final Sample 


411 


416 



Note. 



Effects of sample selection described in 



Table 6. Cepheid Variables: Calculated Quantities 





Period F-Tost 




EiV - I) 


E{B-V) 12 + log(0/H) 


ID 


(days) ratio 


(mag) 


(mag) 


(mag) p/piso (dcx) 



Fl-567 

Fl-773 

Fl-1274 

Fl-1552 

Fl-1619 

Fl-1721 

Fl-1767 

Fl-1890 

Fl-1975 

Fl-2045 

Table to 



40.934 
75.374 
49.414 
42.683 
30.881 
27.260 
21.793 
30.405 
30.419 
22.837 



19.381 
28.559 
82.090 
35.455 
59.668 
44.162 
28.944 
44.678 
45.180 
38.497 



10.350(074) 
10.644(099) 
10.769(101) 
10.626(164) 
10.453(077) 
10.507(072) 
10.303(088) 
10.324(077) 
10.384(072) 
10.610(068) 



be published in its entirety in machine 



0.020(020) 
0.075(074) 
0.066(065) 
0.102(102) 
0.062(062) 
0.017(017) 
0.008(007) 
0.175(175) 
0.163(163) 
0.051(163) 
readable form. 



015(014) 
054(054) 
048(047) 
074(074) 
045(045) 
013(012) 
006(005) 
128(127) 
119(118) 
037(118) 



0.1696 
0.2551 
0.1131 
0.1897 
0.2160 
0.1573 
0.1569 
0.1097 
0.1620 
0.1740 



8.60 
8.52 
8.65 
8.58 
8.56 
8.61 
8.61 
8.65 
8.60 
8.59 



Note. — Errors in the mean colors and the LMC ralative distance moduli are in units of 0.001 mag. 
Table discussed in i|4.2l 
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Table 7. Error Budget of the Cepheid Distance Scale 



Error source 


LMC relative /io 


This Work MlOl 


Goal^ MlOl ^0 


A. Fiducial galaxy 




NGC 4258 


NGC 4258 


Sla. Distance modulus (sys) 




0.12 




Sib. Distance modulus (ran) 




0.09 




SI. Distance modulus (tot) 




0.15 


0.07 


B. Photometric calibration 








S2a. V zeropoint 


0.02 


0.02 




S2b. I zeropoint 


0.02 


0.02 




S2. Photometry (sys) 


0.05 


0.05 




S3. NGC 4258^ Photometry (sys) . . . 




0.05 




Rl. Photometry (ran) 


0.03 


0.03 


0.03 


R2. NGC 4258'^ Photometry (ran) . . 




0.03 


0.03 


C. Extinction corrections 










u.uz 


u.uz 


n 09 
u.uz 


R4. NGC 4258^ Uncertainty in Ry ■ 




0.02 


0.02 


D. Corrections to LMC Metallicity 








S4. Adopted correction 


0.04 


0.04 


0.04 


S5. NGC 4258^ Adopted correction . 




0.04 


0.04 


R5. Zero-point 


0.02 


0.02 


0.02 


R6. NGC 4258^ Zero-point 






0.02 


Rt- Total random 


0.04 


0.05 


0.06 


St- Total systematic 


0.07 


0.18 


0.09 


Combined error (mag) 


0.08 


0.18 


0.11 


Combined error (%) 


4 


8 


5 



Note. — All errors expressed in magnitudes unless otherwise indicated. Table discussed in §5.31 and 



^ Goal erro rs are calculated using the anticipated improvement in the maser distance determination to 
NGC 4258 dHumphrevs et ahlbood l. 



This work uses the 



Macri et al. 



(j2006l ) C epheid NGC 4258 d istance. Goal work errors are based on 



using the same analysis in this paper on the iMacri et al.l (|2006l ) HST/ACS Cepheid observations. 
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Table 8. Recent Distance Determinations to MlOl 





ivietiioa 




iVilUl flQ 


rjivn.^- relative ^/^q 


Outcr/lnncr Field 






mag 


mag 


mag 




This work 


Ccpheids 


18.41±0.10r±0.13s 


29.04±0.05r±0.18s 


10.63 ± 0.02r ± 0.04s 


inner 


This work 


TRGB 




29.05±0.06r±0.12s 




inner 


Saha et al. (2006) 


Cepheids 


18.54'= 


29.16 ±0.04 


10.62 ±0.04 


inner 


Saha et al. f2006) 


Cepheids 


18.54'= 


29.18 ± 0.08 


10.64 ±0.08 


outer 


Freedman et al. (2001) 


Cepheids 


18.50 


29.13 ± 0.11 


10.63 ±0.11 


outer 


Macri et al. f2001") .... 


Cepheids 


18.50 


29.19 ± 0.08 


10.69 ±0.08 


inner 


Macri et al. (2001) .... 


Cepheids 


18.50 


29.53 ± 0.03 


11.03 ±0.03 


outer 


Stetson et al. C 1998) . . 


Cepheids 


18.50 


29.21 ± 0.17 


10.71 ±0.17 


inner 


Kennicutt et al. (19981 


Cepheids 


18.50 


29.21 ± 0.09 


10.71 ±0.09 


inner 


Kcnnicutt et al. (1998) 


Cepheids 


18.50 


29.34 ±0.08 


10.84 ±0.08 


outer 


Kelson ct al. (1996) . . . 


Cepheids 


18.50 


29.34 ±0.17 


10.84 ±0.17 


outer 


Rizzi et al. (2007) 


TRGB 




29.34 ±0.09 






Sakai et al. (2004) .... 


TRGB 




29.42 ± 0.02 






Jurcevic et al. (2000) . 


RSV 


18.50 ±0.10 


29.40 ±0.16 


10.90 ±0.16 




Feldmeier et al. ( 1996) 


PNLF 




29.42 ± 0.15 







Note. — Table discussed in ^531 JJand ^ 
" TRGB: tip of the red giant branch; RSVired supergiant variable stars; PNLF: planetary nebula luminosity function. 

We simply subtract each measurements adopted LMC Distance Modulus from the MlOl Distance Modulus to obtain the LMC- 
relative Distance Modulus for rough comparisons with our measurements. 

'= Zero poin t also assumes a dis tance modulus to Pleiad s is 5.61 mag and relays on the Baade-Becker-Wesselink distances to 32 
Cepheids in lFouque et al.l ||2003| ) and lBarnes et al.l ||2003|) . 



Table 9. Model Wesenheit P-L Fit Parameters 



Model 














DoF 




mag log (day) 


mag log (day) ^ 


mag 


mag 






Model 1 






_n 1+0.1 




22.21 


+0.01 
-0.01 


278.8 


279 


Model 2 






_o ,+0.1 

■^■^-0.1 


0.681°;?° 


22.22 


+0.01 
-0.01 


263.4 


278 


Model 3 


-3.4l 


1.5 

1.6 


_o 1+0.1 
■^■^-0.1 




22.20 


+0.01 
-0.01 


273.8 


278 


Model 4 


-2.8l 


1.6 

1.6 


_o 1+0.1 

■^■^-0.1 


0.651°;?° 


22.21 


+0.01 
-0.01 


260.1 


277 


Restricted Model 1 










22.20 


+0.01 
-0.01 


265.0 


265 


Restricted Model 2 








0.681°;?° 


22.21 


+0.01 
-0.01 


251.3 


264 


Restricted Model 3 


-3.5l 


1.6 

1.7 


_o 1+0.1 
■^■^-0.1 




22.20 


+0.01 
-0.01 


260.7 


264 


Restricted Model 4 


-2.81 


1.7 

1.7 




0.651°;?° 


22.21 


+0.01 
-0.01 


248.8 


263 



Note. — 

Restricted Models have a 40.0 day maximum period cut imposed. 68% uncertainties and confidence 

intervals are presented. Table discussed in 36]and ^ 

Model 1: constant slope and intercept 

Model 2: linearly varying intercept with a constant slope 

Model 3: linearly varying slope with a constant intercept 

Model 4: linearly varying slope and intercept 
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Table 10. Model Wesenheit P-L Fit Parameters as a function of [0/H] 



Model 










mag log (day ) ~ ""^ dex"-*^ mag log(day)~"'^ 


mag dex~^ 


mag 



Model 1 ■■■ -3.1«;i ■-■ ^^--o.oi 

Model 2 ■■■ -3.ll^;} ~0.75tail±0ms 22.22j;;i[;°J ± 0.04^ 

Models 3.8li;7±0.3s -S.lto.} ± 0.2^ ■■■ 22.20j:o;°i 

Model 4 3.ll};*±0.2s -3.lto.} ± 0.2^ -0.72j;[J;2° ± 0.06^ 22.2llg;°j ± 0.04^ 

Restricted Model 1 ■ ■ ■ -3-Oto } ■ ■ ■ '''' '^'^+°-°^ 

,+0:i _n '7C+0.21 



Restricted Model 2 ■ ■ ■ -S-OIg i -0-76_o 22 ± "-O^s 22.21I^ ;5; ± 0.04s 

Restricted Model 3 3.91};^ ± 0.3s -S.iIq.} ± 0.2^ ■ ■ ■ 22.20j;o;°i 

Restricted Model 4 3.ll};g ± 0.2^ -S.Otg i ± 0.2^ -0.72j;[];2^ ± 0.06s 22.2ll[5;{;j ± 0.04s 



Note. — Same as Table [9] but model fits are now a function of metallicity using Eq|9] 68% uncertainties and 
confidence intervals are presented. Table discussed in ^and il9l 



Table 11. Model Wesenheit P-L Comparison 



Models Compared 


Full Sample 




Restricted Sample 




F-Test Confidence 




F-Test Confidence Ax^ 


Model 2 vs. Model 1 


99.99% 


15.5 


99.98% 13.7 


Model 4 vs. Model 2 


93.6% 


3.2 


89.4% 2.5 



Note. — Restricted Sample has a 40.0 day maximum period cut imposed. Table 
discussed in ijejand 39] 



Table 12. Other Objects 





R.A. 


Dec. 


X 


Y 


Magi 


nitudcs 




12 + logiO/H) 


ID 


(J2000) 


(J2000) 


pixels 


pixels 


V 


I 


p/Piso 


(dex) 


Fl-13 


210.86366905 


54.35792401 


2375.8 


2560.5 


20.763 


20.145 


0.1577 


8.61 


Fl-29 


210.84194069 


54.34349859 


3688.0 


2994.3 


21.116 


20.243 


0.1012 


8.66 


Fl-113 


210.88622020 


54.37159190 


1100.4 


2074.4 


21.601 


20.522 


0.2286 


8.54 


Fl-54 


210.83356976 


54.33960030 


4083.5 


3208.2 


21.381 


20.970 


0.0868 


8.67 


Fl-258 


210.85598559 


54.37583152 


1314.9 


3360.8 


22.239 


21.909 


0.1771 


8.59 


Fl-347 


210.81751077 


54.34294154 


4125.0 


3922.7 


22.660 


22.840 


0.0450 


8.71 


F2-24 


210.75884312 


54.34062464 


1344.9 


595.3 


20.549 


19.616 


0.1160 


8.65 


F2-91 


210.76574754 


54.35390318 


350.0 


685.4 


21.265 


20.443 


0.0919 


8.67 


F2-37 


210.70780831 


54.32347179 


3292.9 


2118.5 


20.917 


20.519 


0.2704 


8.51 


F2-48 


210.69712403 


54.35234385 


1533.4 


3312.8 


21.029 


21.079 


0.2601 


8.52 


F2-660 


210.70992763 


54.34032812 


2134.3 


2490.6 


22.718 


21.845 


0.2367 


8.54 


F2-367 


210.70670420 


54.34630827 


1785.8 


2777.3 


22.297 


22.200 


0.2392 


8.53 



Note. — R.A. and Dec. are in units of degrees. Astrometry presented is referenced to the HST headers of 
j9o401010.drz.fits for Field 1 and j9o413010.drz.fits for Field 2. Table discussed in ^ 
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Table 13. Other Possible Periodic Objects 





Period 


R.A. 


Dec 


X 


Y 


Magnitudes 




12 + log (O/H) 


ID 


(days) 


(J2000) 


(J2000) 


pixels 


pixels 


V I 


p/piso 


(dex) 


Fl-574 


0.508 


210.82228318 


54.34545471 


3880.3 


3808.9 


22.806 22.747 


0.0516 


8.70 


F2-271 


10.779 


210.75024732 


54.35297493 


656.0 


1263.2 


21.998 21.962 


0.1291 


8.63 



Note. — R.A. and Dec. are in units of degrees. Astrometry presented is referenced to the HST headers of 
j9o401010.drz.fits for Field 1 and j9o413010.drz.fits for Field 2. Table discussed in ^ 



Table 14. Pixel Coordinate Conventions 



Software Coordinate of Center of Lower Left Pixel 

X pixel y pixel 

DOLPHOT 0.5 0.5 

IRAF 1.0 1.0 

DAOPHOT 1.0 1.0 

ISIS 0.0 0.0 

DS9" 1.0 1.0 



Note. — Table discus sed in ^ 
" Ijove Mande] j^Ooi) 
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Fig. 1. — SDSS image of MlOl showing the two fields observed with ACS/WFC (rectangular 
shape) which are used in this study and the two fields observed with WFPC2 ('stealth fighter' 
shape) for the Hubble Key Project. For both projects, Field 1 is in the East and Field 2 is in the 
West. Figure discussed in ^ 




Fig. 2.— HST ACS/WFC F555W and F814W 2-band color composite image of the ISIS reference 
images for Field 1 in MlOl. The absolute orientation of the field is shown in Fig. [TJ The ACS/WFC 
field of view is 202" on each side, corresponding to a physical size of 6.3 kpc. Figure discussed in 



Fig. 3.— HST ACS/WFC F555W and F814W 2-band color composite image of the ISIS reference 
images for Field 2 in MlOl. The absolute orientation of the field is shown in Fig. [TJ The ACS/WFC 
field of view is 202" on each side, corresponding to a physical size of 6.3 kpc. Figure discussed in 

mm 
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Fig. 4. — Inner region of Field 2 in MlOl from our F555W ISIS reference image (left) and om' 
F555W ISIS variability map (right). Cepheids from the final sample found in ^ are marked by 
teal annuli. The region is ~ 62" on each side, corresponding to a physical size of ~ 1.9 kpc. Figure 
discussed in ^3.11 
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Fig. 5. — Welch-Stetson variability index ly (IWelch &: Stetsonlll993l ) versus F-band magnitude 
for the objects in Field 1. The black solid line marks a variability index of 0. The red dashed 
line marks the ly > 1-5 criteria used to select variables. Objects with a /y > 40 are marked with 
arrows. The blue lines represent the 68% and 90% contours. Red circles and red arrows mark the 
locations of the final Cepheid sample found in f|H Figure discussed in §4.11 
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Fig. 6. — Cepheid light curves of various periods from Field 1. Blue (bottom): V filter liglit 
curves. Red (top): / filter light curves. Solid line shows the best-fit light curve template from 
Yoachim et al.l (j2009l ) with additional fit amplitude parameters. Figure discussed in 
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Id: F2-19702, P= 3.36 
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Fig. 7. — Cepheid light curves of various periods from Field 2. Blue (bottom): V filter liglit 
curves. Red (top): / filter light curves. Solid line shows the best-fit light curve template from 
Yoachim et al.l (j2009l ) with additional fit amplitude parameters. Figure discussed in ^4.21 
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Fig. 8. — Distribution of Ax^/dof as a function of the template model period for objects in Field 
1. Objects with Ax^/dof > 30 are represented by arrows. Figure discussed in §4.21 



- 42 - 




3 4 6 10 20 40 60 

Period (days) 



Fig. 9. — Distribution of periods in Field 1 (solid) and Field 2 (dashed) for the final Cepheid 
sample shown in Table El Figure discussed in §4.41 
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Fig. 10.— ISIS reference image of the HST ACS/WFC F555W data for Field 1 in MlOl. Cepheids 
from the final sample shown in Table [3] are marked by teal open circles. The area of each circle is 
proportional to the logarithm of the Cepheid's period. The shortest and longest period Cepheids 
have periods of 3.2 days and 75.4 days, respectively. Figure discussed in ^4.4[ 
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Fig. 11.— ISIS reference image oftheHSTACS/WFCF555W data for Field 2 in MlOl. Cepheids 
from the final sample shown in Table [3] are marked by teal open circles. The area of each circle is 
proportional to the logarithm of the Cepheid's period. The shortest and longest period Cepheids 
have periods of 3.0 days and 69.4 days, respectively. Figure discussed in ^4.4[ 
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Fig. 12. — I versus (V—I) color-magnitude diagram for Field 1 (left) and Field 2 (right), displaying 
every fifth so urce. The solid blue lines are the predicted position of the unreddened LMC P- 
L relations of lUdalski et al.l (jl999al ) for Cepheids of periods 4-100 days, shifted to the distance 
modulus fit in Fig. [16] and Fig. [T71 The tick marks are every 0.2 log (days) along the P-L relation. 
The arrow indicates a reddening of E{B — V) = 0.2 mag. Red circles represent the final Cepheid 
sample presented in Table El Figure discussed in §4.41 
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Fig. 13. — V, I and Wesenheit P-L relations for F i eld 1 ( left) and Field 2 (right). The blue solid 
lines are the LMC P-L relations of lUdalski et al.l (jl999al ) shifted to the distance modulus fit in 
Fig. [m and Fig. [T71 Random errors for both period and magnitude are shown. Reddening can be 
seen in the V and / P-L relations. Figure discussed in 
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Fig. 14. — V, I and Wesenheit LMC relative distance moduli of the Cepheids as a function of 
period for Field 1 (left) and Field 2 (right). The vertical dashed blue line marks the minimum 
period cut used in each field to refine the final Cepheid sample used to determine the distance 
modulus. 155 Cepheids in Field 1 and 135 Cepheids in Field 2 have periods larger then the period 
cut. The horizontal solid blue line represents the distance modulus fit in Fig. [16] and Fig. [T71 
Random errors for both period and relative distance modulus are shown. Reddening can be seen 
in the V and / P-L relations for Cepheids with periods above the period cut. Figure discussed in 
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Fig. 15. — Median relative Wesenheit distance moduli as a function of the minimum period cut 
applied. Filled circles represent Field 1. Red crosses represent Field 2 and are shifted to higher 
periods by 0.01 log(days). The solid line and the dashed red line represents the median distance 
modulus for Field 1. Final choices for period cut are shown by open circles. 155 Cepheids in Field 
1 and 135 Cepheids in Field 2 have periods larger then the period cut. The dashed line represents 
the median distance modulus at the period cut for Fields 1 and 2, respectively. Figure discussed 
in 
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Fig. 16. — Relative Wesenheit distance moduli as a function of period for Field 1. Random 
errors are shown for Cepheids with periods larger then the period cut. Figure discussed in §5.11 
Filled Circles: 155 Cepheids with periods larger then the period cut. Solid line: Median 
relative distance modulus. Dotted line: 1 a dispersion about the median in final Cepheid sample. 
Dashed line: Represents period cut. Red Symbols: Cepheids removed from various cuts. Stars: 
Amplitude ratio cut. Triangles: Blue blends cut. Crosses: Red blend or large extinction cut. 
Squares: Population II Cepheid or egregious outlier cut. Red Circle: Period cut. Circled Red 
Point: Iterative sigma clipping. 
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Fig. 17. — Relative Wesenheit distance moduli as a function of period for Field 2. Random 
errors are shown for Cepheids with periods larger then the period cut. Figure discussed in §5.11 
Filled Circles: 135 Cepheids with periods larger then the period cut. Solid line: Median 
relative distance modulus. Dotted line: 1 a dispersion about the median in final Cepheid sample. 
Dashed line: Represents period cut. Red Symbols: Cepheids removed from various cuts. Stars: 
Amplitude ratio cut. Triangles: Blue blends cut. Crosses: Red blend or large extinction cut. 
Squares: Population II Cepheid or egregious outlier cut. Red Circle: Period cut. Circled Red 
Point: Iterative sigma clipping. 
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Fig. 18. — E{V — I) as a function of period for Field 1. Filled Circles: 155 Cepheids with periods 
larger then the period cut. Solid line: E{V — I) values where the large extinction and blue blend 
cuts were apphed. Vertical Dashed line: Represents period cut off. Red Symbols: Cepheids 
removed from various cuts. Stars: Amplitude ratio cut. Triangles: Blue blends cut. Crosses: 
Red blend or large extinction cut. Squares: Population 11 Cepheid or egregious outlier cut. Red 
Circle: Period cut. Circled Red Point: Iterative sigma clipping. Figure discussed in ^5.11 
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Fig. 19. — E{V — /) as a function period for Field 2. Filled Circles: 135 Cepheids with periods 
larger then the period cut. Solid line: E{V — I) values where the large extinction and blue blend 
cuts were apphed. Vertical Dashed line: Represents period cut off. Red Symbols: Cepheids 
removed from various cuts. Stars: Amplitude ratio cut. Triangles: Blue blends cut. Crosses: 
Red blend or large extinction cut. Squares: Population 11 Cepheid or egregious outlier cut. Red 
Circle: Period cut. Circled Red Point: Iterative sigma clipping. Figure discussed in ^5.11 
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Fig. 20. — Difference between median and mean relative Wesenheit distance moduli as a function 
minimum period cut applied. Top: Field 1. Bottom: Field 2. Solid red line is at to guide the 
eye. Errors estimates shown are the quadrature addition of the scatter about the median and the 
scatter about the mean. Figure discussed in ^5.1[ 
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Fig. 21. — Metalicity dependence of the relative Wesenheit distance modulus. Filled circles and 
blue crosses represent Cepheids in the refined sample in Fields 1 and 2, respectively. The solid red 
best-fit line is show with a slope of —0.84 ± 0.22,. it 0.07^ mag dex~^. The dotted red best-fit 
constant is shown for reference. The dashed red best-fit line to a restricted sample of Cepheid with 
12 + log(0/H) < 8.65 dex is shown with a slope of -0.66 ± 0.25^ ± 0.05^ mag dex"!. At 12 + 
log(0/H) = 8.5 dex, the LMC relative distance modulus of MlOl is A/iL^c = 29.04 ± 0.02^ ±0.04^. 
The error bar on the open circle in lower right displays the average individual random uncertainties. 
Figure discussed in §5.21 
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Fig. 22. — Residuals about the best-fit metallicity relation shown in Fig. [2TJ Filled circles and 
blue crosses represent Cepheids in the refined sample in Fields 1 and 2, respectively. Dashed red 
line shown at Afij^^Q [O — C] = to aid the eye. Figure discussed in 
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Fig. 23. — Wesenheit P-L relations with equal numbers of Cepheids in bins of deprojected radius. 
Filled circles represent Cepheids selected in i j4.4l with an additional minimum period cut of 14 
days. The dotted black line indicates fit Wesenheit P-L relation with no variation as function 
of deprojected radius (Model 1). The dashed red line represents the best-fit Wesenheit P-L with 
constant a slope and a linearly varying zero-point as a function of deprojected radius (Model 2). 
The solid blue line is the best-fit Wesenheit P-L allowing linear variations in the slope and zero- 
point as a function of deprojected radius (Model 4). Cepheids were removed by iterative sigma 
clipping from Model 1 and are represented by red points encircled by red rings. The mean Cepheid 
deprojected isophotal radius fraction is shown for each frame. Fit parameters shown in Table [H 
Figure discussed in f|6] and ^ 
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Fig. 24. — The best-fit slope of the Wesenheit P-L relation as a function of fractional isophotal 
radius. The dotted black line indicates fit Wesenheit P-L relation with constant slope and zero- 
point (Model 1). The dashed red line represents fit Wesenheit P-L with constant slope and linearly 
varying zero-point as a function of deprojected radius (Model 2). The dot dashed purple line 
indicates fit Wesenheit P-L with constant zero-point and linear varying slope as a function of 
deprojected radius (Model 3). The dashed red line represents fit Wesenheit P-L with linear varying 
slope and zero-point as a function of deprojected radius (Model 4). The filled black dots represent 
the slopes of the individually fit Wesenheit P-L relations for smaller samples shown in Fig. [23] with 
error bars determined using bootstrap re-sampling. It is important to note that these individually 
fit sub-samples are plotted for comparison purposes only. Figure discussed in ^ 




Fig. 25. — Same as Fig. [23] but showing the Wesenheit P-L zero-point as a function of fractional 
isophotal radius. Figure discussed in ^ 
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Fig. 26. — Same as Fig. [23] but for the sample Cepheids with a maximum period cut of P < 40.0 
days. Fit parameters shown in Table [9l Figure discussed in ^}6] and ^ 
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Fig. 27.— Determin ation of the T TRGB magnitude for MlOl where T = I- 0.20[{V -I)- 1.5] 



(jMadore et al.M2009l ). Left: T vs. {V — I) color-magnitude diagram for sources 4.75 from the 



galactic center of MlOl. Right: Edge detection function E(m)y'(j){m), where the maximum value 
at 25.00. Solid red line: Represents the maximum value of the edge detection function. Dashed 
red lines: Represent estimated uncertainty of the edge detection based on bootstrap re-sample 
described in fJT) 
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Fig. 28. — Light curves of interesting luminous variable objects. Blue (top): V filter light curves. 
Red (bottom): I filter light curves. Photometric errors are shown but usually smaller then the 
points. Figure discussed in SjSJ 
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Fig. 29. — Object Fl-574 (top) and F2-271 (bottom) which show possible periodic behavior with 
P = 0.508 days and P = 10.78 days respectively. Absolute magnitudes on the right side of the 
figure assume our Cepheid distance modulus and make no reddening correction. Figure discussed 
in 381 
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Fig. 30. — / vs. {V — I) color-magnitude diagram for MlOl for Field 1 (left) and Field 2 (right) 
displaying every fifth source. Red circles represent luminous variable object seen in Fig. [28] and 
Fig. [29j Red annuli mark the possible periodic objects. MS labels the location of the main sequence. 
BSG labels the location of blue super giants. RSG labels the location of red super giants. RG labels 
the location of red giants. Absolute magnitudes on the right side of the figure assume our Cepheid 
distance modulus and make no reddening correction. Figure discussed in ^ 



